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1 General 


1.1 Background 

Upper Trishuli-1 HEP is constructed to utilize the hydro potential of Trishuli River near 
Dunche and Haku. The construction site is situated in a rugged mountain between 
Tibet in China and the Himalaya Mountain, about 70 km away from capital of Nepal, 
Kathmandu, towards the northwest. The project is a run-of-river development scheme. 
The headwork site is at about 280m downstream from the confluence of Bhotekoshi 
and Trishuli rivers. The powerhouse site is 500m upstream from Mailun Dobhan 
where Haku VDC is placed. The installed capacity of Upper Trishuli-1 HEP is 216 MW. It 
is a year-round power plant capable of generating about 1,533.4 GWh electricity 
yearly on average. 

The Adit No.3 is located on the left bank of headrace tunnel, and the headrace tunnel 
chainage at intersection area is H+6562.747m. The cross section of headrace tunnel is 
horseshoe shape and 4.2mx5.6m size. The elevation of the inlet is EL. 1194.28m, and 
the elevation of outlet is EL.1212.36m. The length of Adit No.3 is 309.575m and 
longitudinal slope is 5.86%. 

According to the Contractor’s interpretation and preliminary geological survey data, 
along the tunnel axis, there are interbedded mica schist and quartz schist are 
distributed and they are classified as Type IV,III, I] as per Q-method from inlet to 
outlet, the buried depth from inlet to outlet is about 50m to 300m, averaged 200m 
approximately. It is speculated that the whole tunnel is mostly dry, individual joints 
surface is wet. No potential continuous sliding surface is found on the tunnel slope. 
The longitudinal profile of Adit No.3 is shown in Figure 1.1. The preliminary excavation 


cross section and preliminary support pattern are shown in Figure 1.2~1.7. 
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Figure 1.2 Typical section of adit No.3 (rock class |) 


Page 2 of 50 


DOOSAN 


UPPER TRISHULI-1 HEP (216MW) 


Ch EB 5 


POWERCHINA 


€ 


Surrounding Rock Stability Calculation of Adit No.3 


Rockdowel 
D25mm, L=3m 
(where necessary) 


Weep hole 


L=0.5m,045mm 
(where necessary) 


Plastic fiber reinforced shotcrete —__ 


thickness 50mm 


4200 


200, || 50 


50 | 1200 


Figure 1.3 Typical section of adit No.3 (rock class II) 
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Figure 1.4 Typical section of adit No.3 (rock class III) 
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Figure 1.5 Typical section of adit No.3 (rock class IV) 
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Figure 1.6 Typical section of turning adit (rock class III) 
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Figure 1.7 Typical section of Junction (rock class II) 


1.2 Design Objective 

This calculation only includes rock class I~IV, and the rock class V shall be considered 
separately according to the actual situation. 

This calculation is provided for detailed design stage and aim to confirm the objectives 
which are as following: 

(1) to generally confirm the stability of the tunnel during construction; 

(2) to confirm (or adjust) the previously suggested support measures in the 


construction drawings. 
1.3 Purpose and Contents 


(1) Contract documents, Employer’s Requirement; 

(2) Tunnels and shafts in Rock (EM 1110-2-2901); 

(3) Rock Foundations (EM 1110-1-2908); 

(4) Building Code Requirements for Structural Concrete and Commentary (ACI 318) 


(5) Vlachopoulos, N. and Diederichs, M.S. (2009). Improved longitudinal displacement 
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profiles for convergence-confinement analysis of deep tunnels. Rock Mechanics and 
Rock Engineering, 42(2), 131-146; 

(6) The 2008 Kersten Lecture Integration of geotechnical and structural design in 
tunneling; To be presented at the opening keynote address, by Dr Evert Hoek, at the 
University of Minnesota 56 Annual Geotechnical Engineering Conference to be held 
in Minneapolis on 29 February 2008; 

(7) Handbook Using the Q-system, Rock mass classification and support design; 

(8) Paper: Rock-Support Interaction analysis for tunnels in weak rock masses; 

(9) https://www.rocscience.com/help/rs2/turorials/rs2 _tunnel lining design.htm; 

(10) Drawing: Excavation and initial support drawing of adit No.3 (Drawing No. 
UT1-C-150-CVL-DG-43004). 


2 Theoretical Background of Calculations 


2.1 Software 


Software Phase 2D will be applied in the calculation. Phase2 is a powerful 2D finite 
element program for soil and rock applications. Phase2 can be used for a wide range 
of engineering projects including excavation design, slope stability, groundwater 
seepage, probabilistic analysis, consolidation, and dynamic analysis capabilities. 

Software UNWEDGE is an interactive software for the stability analysis of 
three-dimensional wedge formed by structural discontinuity and underground 
excavation, which is used to analyze the underground excavation problem with 
discontinuous structural plane in rock mass. UNWEDGE calculates the safety factor of 
potentially unstable wedge, and can analyze the influence of support system on 


wedge stability. 
2.2 Modelling Conditions, Material Law, Constitutive Model 


The modelling conditions are plane strain. The material law is elasto-plastic combined 
with strain-softening (for rock). The rock mass is modelled as a continuum, 


homogeneous and isotropic. The failure criteria is the Mohr-Coulomb failure criteria. 
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What needs to be emphasized is that such continuum model with homogeneous and 
isotropic material is a model, only, and does not necessarily reflect reality. In reality, 
the displacements occur in blocky rock masses always along discontinuities (joints, 
foliation etc.), which leads to loosening of the rock mass. The determined plastic zone 
in the continuums model represents an area of the rock mass where rock blocks were 
displaced to each other (= loosening of rock mass). In reality, it does not represent an 


area where the material was subject to plastic flow, as in soil. 


2.3 Calculation Methodology 


Convergence-confinement analysis for tunnelling is a standard approach for 
preliminary analysis of anticipated wall deformation and support design in squeezing 
ground. 

Convergence-confinement analysis is a widely used tool for preliminary assessment of 
squeezing potential and support requirements for circular tunnels in variety of 
geological conditions and stress states. An analytical plasticity solution such as that 
developed by Carranza-Torres and Fairhurst (2000) is applied to a circular opening in 
an isotropic stress field. An internal pressure, initially equal to the in-situ stress is 
applied on the inside of the excavation boundary. The pressure is incrementally 
relaxed until the excavation boundary condition is that of zero normal stress. The 
extent of plastic yielding and thereby, the boundary deformation is calculated at each 
stage of the process. The result is a continuous representation of the 
deformation-internal pressure relationship for the tunnel given a particular material 
strength, deformability, dilation and stress state. 

The internal press is, of course not a representation of reality but rather a surrogate 
for the effect of the gradual reduction of the radial resistance provided by the initially 
present tunnel core (material inside the tunnel boundary) transitioning to an exposed 


boundary with a progressively distant tunnel face and ultimately a long open tunnel 
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with plane strain conditions. The internal pressure that is coupled with a given 
boundary displacement is a measure of the amount of support resistance required to 
prevent further displacement at that point in progressive tunnelling model. 

Elastic deformation of the rock mass starts about two diameters ahead of the 
advancing face and reaches its maximum value at about two diameters behind the 
face. At the face position about one third of the total radial closure of the tunnel has 
already occurred and the tunnel face deforms inwards as illustrated in Fig. 2.1. 
Whether or not these deformations induce stability problems in the tunnel depends 
upon the ratio of rock mass strength to the in-site stress level. The tunnel 
displacement profiles for the roof and the invert of a tunnel are shown schematically 


in the Fig. 2.2. 
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Figure 2.1 The shape of the deformed tunnel profile and displacement vectors 
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Longitudinal Displacement Profile __ 


Radial displacement 
reaches its final value 
at about two tunnel 
diameters behind the 
face 


Inward deformation 
of tunnel face 


— Radial displacement reaches 
about one third of its final 
value at the tunnel face 


S 


Direction of 
face advance 


Displacements start about two 
tunnel diameters ahead of the face 


Figure 2.2 pattern of elastic deformation in the rock mass surrounding an advancing tunnel 
The equations can be found in the Kersten Lecture, appendix 1(Hoek et.al.,2008) 
According to the formula listed in geological rock classification and documentation 
“Integration of geotechnical and structural design in tunneling (Hoek, E., 
Carranza-Torres, C., Diederichs, M.S. and Corkum, 2008)”. The calculation formula of 


supporting time is as follow: 


u, 1 ase 
Uax 3 
P=R /R 
P t 
oX 
a 
t 
Uod 
1g" when¥ <0 
u M max 
7 u 3d. 
= 1-(1- ——)e 2R when¥ > 0 
max 
where, 
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P —The maximum plastic zone depth, m; 


u e 
° —The radial displacement of tunnel face, m; 


R, 


—Tunnel radius, m; 


X —Distance from the initial support to the tunnel face, m; 


Hirax —The maximum radial displacement, m. 

With these analytical equations it is possible to estimate the amount of displacements 
at the specified location of support installation, if the plastic radius and displacement 
far from the tunnel face are known. Therefore, to get the displacements and the 
radius of the plastic zone far away from the tunnel face (where the internal pressure is 
zero). It is necessary to calculate numerically with Phase 2 the tunnel without support, 
first. 

With the determined displacement profile along the tunnel and a given location of 
support installation, the modeler selects the corresponding deformation. 

In the next step the tunnel advance is modelled in Phase 2 by reducing the internal 
pressure in the future tunnel until the corresponding deformation is reached, when 
support is installed. The modeler installs now the support in the model. The internal 
pressure is then reduced to zero by factor the field stress zero. Further numerical 


calculation leads to loading of the support. 
2.4 Maximum Unsupported Span 


The maximum unsupported span is calculated by the following formula as per EM 
2901: X=2 ESR Q^0.4 
Where, X is maximum unsupported span, ESR is excavation support ratio, taken as 1.3. 


Q values and calculated and applied maximum unsupported spans are shown in Table 


fa TR DOM HIRAI ERE ICE. -1 


Table iR CRS HY HB Ee FE TL SC. -1 Maximum unsupported span of tunnel excavation 
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Calculated maximum Maximum unsupported 
Rock class Q Value 
unsupported span(m) span for application, X(m) 
| 40 11.4 11 
Il 25 9.4 9 
II 7 5.7 5 
IV 2.5 3.8 3 


2.5 Ground Characteristic Curve 
Define the relationship between stress release ratio Aj and the tunnel deformation ui 


for an advancing tunnel in stress field. A plot of ui versus Aj, the curve is based on the 


assumption that the rock at the tunnel face provides an initial support pressure equal 
to the in-situ stress. As the tunnel face advances and the face moves away from the 
section under consideration, the support pressure gradually decreases until it reaches 
zero at some distance behind the face. The characteristic curve of class I~IV are shown 


as below. 
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Figure 2.3 Characteristic curve of class | 
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Figure 2.4 Characteristic curve of class II 
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Characteristic curve of Class III 


Total Displacement ui [m 
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Figure 2.5 Characteristic curve of class III 


Characteristic curve of Class IV 


Total Displacement ui [m] 


0% 20% 40% 60% 80% 100% 
Stress Release Ratio Ai 


Figure 2.6 Characteristic curve of class IV 
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2.6 Longitudinal displacement profile (LDP) for tunnel 


The Longitudinal Displacement Profile (LDP) is required in order to establish the 


relative position of the tunnel face and the section under consideration. 


According to the LDP equations, presented above, the LDP curve of every rock mass 


class are shown in the table below. 


The displacement of vertex (mm) 
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| 
2a 
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LDP Curve of Class | 
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Figure 2.7 Longitudinal Displacement Profile curve of class | 
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Figure 2.8 Longitudinal Displacement Profile curve of class II 


LDP Curve of Class III 
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Figure 2.9 Longitudinal Displacement Profile curve of class III 


LDP Curve of Class IV 
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Figure 2.10 Longitudinal Displacement Profile curve of class IV 


According to Characteristic Curve and LDP Curve of class I~IV, it is determined to adopt the 


relaxation of stress at the supporting time shown as below. 
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Table fix! CRY AKA TRE PESKY MH. -2 Relaxation of rock mass 


Supporting Plastic zone | Vertex u 
Stress Relaxation 
Rock Class time X radius Rp displacement Umax 
Closure ratio | Percentage (%) 

(m) (m) Umax (mm) 
l 11 3.62 0.5 0.99 95 
II 9 4.44 1.5 0.97 95 
III 5 5.15 9.1 0.83 90 
IV 3 8.84 62.5 0.53 85 


Relaxation refer to that (100-Relaxation)% of surrounding rock stress is undertaken by rock and 


support, Relaxation % release during excavation. The stress relaxation percentage has already 


considered the margin of safety. 


3 Properties of Materials 


3.1 Rock Mass Properties 


The mechanical parameters of surrounding rock are shown in Table 3.3. 


Table 3.1 Mechanical parameter of surrounding rock 


Uniaxial Rock mass/rock mass Uniaxial 
Elastic 
Class of compressive Poisson’s shear (shearing) compressive | Tensile 
Unit modulus 
surrounding strength ratio of strength strength of |Strength of 
weight Of rock 
rock (UCS) rock mass rock mass, | rock mass 
mass Peak Residual 
of intact rock (UCS) 
Y Oc Em u O peak Cpeak Ores Cres Ocrm Otm 
(kN/m?) (MPa) (GPa) (°) (MPa) (°) (MPa) (MPa) (MPa) 
I 27.0 60 30.0 0.20 50 5.0 | 40 | 3.33 27:5 2.75 
O 27.0 60 20.5 0.25 48 2.5 38 | 1.67 13.0 1.3 
26.5 60 he 0.27 45 1 36 | 0.66 4.2 0.42 
IV 26.0 60 2.3 0.30 35 0.5 28 | 0.33 1:9 0.19 


Note: 1) Residual strength parameters shall be reduced by 1/3 for the cohesion and 20% for the 


friction angle; 2) Dilation angle was estimated at 1/3 of the friction angle; 3) Tensile strength of 


rock mass is taken as 10% of uniaxial compressive strength of rock mass. 
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3.2 Shotcrete 


a) Plastic fiber shotcrete: unit weight: 24.0kN/m3; 
b) The Elastic moduli are calculated as formula E-=4700*f,’*0.5 (ACI 318, section 


19.2.2 Modulus of elasticity). 


Table 3.2 Shotcrete strength and elastic modulus 


Compressive strength f (MPa) Elastic modulus E.(MPa) Poisson ratio 
25 23500 0.2 
30 25743 0.2 


3.3 Rock Dowels 


Rock dowel type: Cement mortar full length bonded dowel 

1) Yield strength: SOOMPa 

2) Compressive/Tensile strength: 545MPa 

3) Elastic modulus: 210,000MPa 

4) Tensile strength of rock bolt steel: SOOMPa/1.5=333MPa, 

And the value 1.5 is safety factor. 
Ultimate bearing capacity of rock dowel is calculated from formula below: 
F=0.5rtdLtuit 

Where, 

d =Effective diameter of borehole; 

L = Length of grouted portion of anchor; 

t =Working bond strength; 

Tut= The ultimate bond strength at failure, the ultimate bond stress is often taken as 
1/10 of the uniaxial compressive strength of the rock or grout (whichever is less) 
(Littlejohn 1977) up to a maximum value of 4.2 MPa. 


Table 3.3 The ultimate bearing capacity calculation of rock dowel 


Rock classification Unit I IIl IV 
UCS MPa 60 60 60 
Rock dowel diameter mm 25 25 25 
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Rock dowel area mm? 491 491 491 
Ultimate bond strength of rock dowel Tutt kPa 4200 4200 4200 
Effective diameter of rock dowel borehole d | mm 50 50 50 
; 0.5x50x4200/10%= 
Bond strength of unit length F=0.5rtdTutt MN/m 0.33 aise 0.33 
Tensile strength of rock dowel steel MPa 333 333 333 
f 491x333/10°= 
Tensile peak value of rock dowel MN 0.164 aes 0.164 


3.4 Lattice Girder a 


nd Reinforcing Bar 


The mechanical parameters of reinforcing bar are shown as below. 


Table 3.4 Mechanical parameter for reinforcing bar 


Yield strength (MPa) 


Tension strength(MPa) 


Elastic modulus E.(MPa) 


Poisson ratio 


500 


545 


210,000 


0.3 


3.5 Steel support 


MB150 H-beam type steel rib is used. And the mechanical parameters of steel support 


are shown as below. 


Table 3.5 Mechanical parameter for steel support 


Yield strength (MPa) 


Tension strength(MPa) 


Elastic modulus E.(MPa) 


Poisson ratio 


250 


400 


210,000 


0.3 


4 Initial Proposed Tunnel Supports 


4.1 Tunnel Support 


of Excavation Class 


The initial proposed supports for each rock type are shown as below. 


Table Hix CRS HLA HSE RESIN CF. -3 Initial proposed supports for each rock type 


Initial proposed supports 
Item Unit 

l II III IV 
Shotcrete thickness mm Spot 50 80 100 
Shotcrete strength MPa 21 21 21 21 
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Rock dowel length m 3.0 3.0 3.0 3.0 
Rock dowel spacing within each row m spot spot 2.0 1.5 
Rock dowel row spacing m spot spot 2.0 1.5 
Steel support/Lattice girder spacing m spot 
4.2 Safety Factor 


(1) The support safety factor for shotcrete and rock dowel are 1.5. 


(2) The minimum safety factor for wedges analysis is 2.0. 


5 Model Set-up, Boundary and Loading Conditions, Stress Initiation 


For simplicity, the model will not extend to the ground surface, and the boundary of 
model is set 10 times tunnel width/height, but to apply corresponding stresses at the 
model’s upper boundary. 

For the case of gravitational loading, only, boundary and loading conditions to be set 
as shown below. This applies for tunnels in the UT1-Project that are oriented 
NNE-SSW (+/- 45°), as NE-SSW is generally orientation of the maximum horizontal 
stress Op in the Upper Himalaya, which is in the plane-strain Phase2-model the 
(“irrelevant”) out-of-plane stress. Hence, initial stress conditions in the model can be 
calculated by gravitational loading, or by specifying the relative ground surface 


elevation. 


The finite element model for rock class |, Il, II, IV is shown as SiR! ARFRBI S| FARR. , 


and the corresponding gravitational loading is shown in the i! AA 5| AR. . 
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| | apply | | | | YH | ] 
stress i i j 
yxh f ) 


x fixed, x fixed, i 
y free y free ; 


| t l 
| ; | 


x 


x fixed, 
y fixed 


; 
i 


Figure 5.1 Finite element model and gravitational loading 


Table 5.1 gravitational loading for surrounding rock 


Tunnel Rock Overburden 
Rock class/ Ss os 
, buried unit gravitational | Application location of Max. buried depth 
Location ; ; 
depth weight loading 
m kN/m? MPa 
Refer to rock class II location 
I 300 27 8.1 
Whole tunnel 
I 300 27 8.1 
. Class Il 
Junction 300 27 8.1 
Whole tunnel 
I 300 26.5 7.95 
. Class Ill 
Turning bay 150 26.5 2.2 
200 26 5.20 Refer to rock class III location 
IV 150 26 3.60 Intermediate value 
100 26 2.60 Inferred location geologically 


6 Tunnel Modelling and Results 


6.1 Rock Class | (MAX. Buried Depth 300m) 


For surrounding rock class | without support, adopting the parameters listed in Section 


3, and the calculated displacement and plastic zone are shown as below. 
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Total 
Displacement 
m 

0.00e+000 


1.40e-004 
2.80e-004 
4.20e-004 
5.60e-004 
7.00e-004 
8.40e-004 
9.80e-004 
1.12e-003 
1.26e-003 
1.40e-003 
1.54e-003 
1.68e-003 


5.10e-004 


1.68e-003 


Figure 6.1 Plastic zone and total displacement without tunnel support 


It can be seen that the maximum displacement is 0.5mm in roof and 0.9mm in invert, 


and almost no plastic zone. In terms of stability, the plastic fiber shotcrete may be 


required when necessary. 


6.2 Rock Class II (MAX. Buried Depth 300m) 


For rock class Il, according to LDP curve and Characteristic curve presented in section 


2, it is determined to choose 9m as the max. unsupported distance at support 


installation. The displacement without tunnel support, the displacement after support 


installation, plastic zone and axial force of rock dowels, the safety factor of plastic fiber 


shotcrete/ steel rib (including shear and moment) are shown as following. 


Total 
Displacement 
m 

0.00e+000 


Figure 6.2 Total displacement without tunnel support 
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Total 
Displacement 


m 


Yielded 
Elements 
[per-cent] 


0.00e+000 1.50e-003 


3.00e-004 
6.00e-004 
9.00e-004 
1.20e-003 
1.50e-003 
1.80e-003 
2.10e-003 


2.40e-003 
2.70e-003 
3.00e-003 
3.30e-003 
3.60e-003 


2.25e-003 


Figure 6.3 Total displacement after support installation 
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Figure 6.4 plastic zone and axial force of rock dowels after support installation 
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YI 


= Stagell 
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Support Element: plastic fiber shotcrete 


Moment M- MNM 


Figure 6.5 the safety factor of plastic fiber shotcrete 
6.3 Rock Class III (MAX. Buried Depth 300m) 


For rock class Ill, according to LDP curve and Characteristic curve presented in section 


2, it is determined to choose 5m as the max. unsupported distance at support 


installation. The displacement without tunnel support, the displacement after support 


installation, plastic zone and axial force of rock dowels, the safety factor of plastic fiber 


shotcrete/ steel rib (including shear and moment) are shown as following. 


Total 
Displacement 
m 


0.00e+000 


9.06e-003} 


Figure 6.6 Total displacement without tunnel support 
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Total 
Displacement 
m 

0.00e+000 


9.00e-004 

-80e-003 
2.70e-003 
3.60e-003 
4.50e-003 
5.40e-003 
6.30e-003 
7.20e-003 
8.10e-003 
9.00e-003 
9.90e-003 
1.08e-002 


Yielded 
Elements 
[per-cent] 
0 


50 


100 


5.85e-003 


Figure 6.7 Total displacement after support installation 


SSS 


0.159MN Axial Force 


Figure 6.8 plastic zone and axial force of rock dowels after support installation 
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Support Element: plastic fiber shoterete 


Figure 6.9 the safety factor of plastic fiber shotcrete 


6.4 Rock Class IV (MAX. Buried Depth 100m) 
For rock class IV, According to LDP curve and Characteristic curve presented as 
following, it is determined to choose 3m as the max. unsupported distance at support 


installation. 


Characteristic curve of Class IV 
1.60E-02 


1.40E-02 


1.20E-02 


1.00E-02 


8.00E-03 


6.00E-03 


Total Displacement ui [m] 


4.00E-03 


2.00E-03 


0.00E+00 


0% 20% 40% 60% 80% 100% 
Stress Release Ratio Ai 


Figure 6.10 Characteristic curve 
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14.005 


12.-00— 


10-00 


The displacement of vertex (mm) 


. 0-00 
-10.0 -5.0 0. 0 5. 0 10.0 15.0 


The distance from tunnel advancing face (m) 


| 
ps 
=) 
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Figure 6.11 Longitudinal Displacement Profile curve 


The displacement without tunnel support, the displacement after support installation, 
plastic zone and axial force of rock dowels, the safety factor of plastic fiber shotcrete/ 


steel rib (including shear and moment) are shown as following. 


Total 
Displacement 
m 

0.00e+000 


1.40e-003 
2.80e-003 
4.20e-003 
5.60e-003 
7.00e-003 
8.40e-003 


9.80e-003 
1.12e-002 
1.26e-002 
1.40e-002 
1.54e-002 


1.68e-002 


1.68e-002 


Figure 6.12 Total displacement without tunnel support 
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“{7.80e-003} 


Total 
Displacement 
m 

0.00e+000 


1.30e-003 
2.60e-003 
3.90e-003 
5.20e-003 
6.50e-003 
7.80e-003 
9.10e-003 
1.04e-002 
1.17e-002 
1.30e-002 
1.43e-002 


1.56e-002 


Figure 6.13 Total displacement after support installation 
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[per-cent] 
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Figure 6.14 plastic zone and axial force of rock dowels after support installation 
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Concrete: 0.1m [Togni 


Thrust N- MN 


Support Element: plastic fiber shotcrete 


Figure 6.15 the safety factor of plastic fiber shotcrete 


6.5 Rock Class IV (MAX. Buried Depth 150m) 


For rock class IV, According to LDP curve and Characteristic curve presented as 
following, it is determined to choose 3m as the max. unsupported distance at support 


installation. 


Characteristic curve of Class IV 
4.00E-02 


3.50E-02 


3.00E-02 


2.50E-02 


2.00E-02 


1.50E-02 


Total Displacement ui [m] 


1.00E-02 


5.00E-03 


0.00E+00 
0% 20% 40% 60% 80% 100% 


Stress Release Ratio Ai 


Figure 6.16 Characteristic curve 


Page 28 of 50 


DOOSAWN | UPPER TRISHULI-1 HEP (216MW) € tHg 


POWERCHINA 


Surrounding Rock Stability Calculation of Adit No.3 


LDP Curve of Class IV 


The displacement of vertex (mm) 


-15. 0 -10. 0 -5. 0 00 5.0 10.0 15.0 
The distance from tunnel advancing face (m) 
Figure 6.17 Longitudinal Displacement Profile curve 
Some support system is considered: 
(1) The shotcrete thickness is increased to 0.12m; 
(2) The shotcrete strength is increased to 30MPa; 
(3) The spacing of Steel support is taken as 1m; 
(4) The spacing of rock dowel is decreased to 1m; 
The displacement without tunnel support, the displacement after support installation, 
plastic zone and axial force of rock dowels, the safety factor of plastic fiber shotcrete/ 


steel rib (including shear and moment) are shown as following. 
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Total 3.37e-002 


Displacement 


0.00e+000 
4.00e-003 
8.00e-003 
1.20e-002 
1.60e-002 
2.00e-002 
2.40e-002 
2.80e-002 
3.20e-002 
3.60e-002 
4.00e-002 
4.40e-002 
4.80e-002 


4.20e-002| 


Figure 6.18 Total displacement without tunnel support 


Total 
Displacement 
m 

0.00e+000 


3.00e-003 
6.00e-003 
9.00e-003 
1.20e-002 
1.50e-002 


1.80e-002 a } 1 $ 
2.10e-002 | | 
2.40e-002 i 

- 1.95e-002 
2.70e-002 2.10e-002 [eseo] 
3.00e-002 D n 


3.30e-002 


3.60e-002 3.30e-002 


Figure 6.19 Total displacement after support installation 
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Yielded 

Elements 

[per-cent] 
0 


50 


Figure 6.20 plastic zone and axial force of rock dowels after support installation 


Concrete: 012m Fag 
. Sap 


Support Element: plastic fiber shoterete 


Figure 6.21 the safety factor of plastic fiber shotcrete 
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6.6 Rock Class IV (MAX. Buried Depth 200m) 


For rock class IV, according to LDP curve and Characteristic curve presented in section 


2, it is determined to choose 3m as the max. unsupported distance at support 


installation. 


Some support system is considered: 


(1) The shotcrete thickness is increased to 0.16m; 


(2) The shotcrete strength is increased to 30MPa; 


(3) The spacing of Steel support is taken as 1m; 


(4) The spacing of rock dowel is decreased to 1m; 


(5) The length of rock dowel is increased to 3.5m. 


The displacement without tunnel support, the displacement after support installation, 


plastic zone and axial force of rock dowels, the safety factor of plastic fiber shotcrete/ 


steel rib (including shear and moment) are shown as following. 


Total 
Displacement 
m 

0.00e+000 


7.00e-003 
1.40e-002 
2.10e-002 
2.80e-002 
3.50e-002 
4.20e-002 
4.90e-002 
5.60e-002 
6.30e-002 
7.00e-002 
7.70e-002 
8.40e-002 


6.41e-002 


5.95e-002 


7.70e-002 


Figure 6.22 Total displacement without tunnel support 


Page 32 of 50 


DOOSAWN | UPPER TRISHULI-1 HEP (216MW) q ChE) Fa 5 


POWERCHINA 


Surrounding Rock Stability Calculation of Adit No.3 


Total 
Displacement 

m 

00e+000 


00e-003 
20e-002 
80e-002 
40e-002 
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60e-002 , 
20e-002 i ~ 
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6.30e-002 


Figure 6.23 Total displacement after support installation 


Page 33 of 50 


POWERCHINA 


DOOSAWN | UPPER TRISHULI-1 HEP (216MW) € ChE) Fa 5 


Surrounding Rock Stability Calculation of Adit No.3 


Yielded 7 ; 
" NNN 
Elements u LRNVZMNINial Force] 
[per-cent] (A> S x Sy 
0 w > aS 
[à 
my 
= 
A 
m: 
= 
=} 
E 
= 
= 
50 = 
= 
E 
= 
mi 
Ll 
E 
=| 
me 
100 


Figure 6.24 plastic zone and axial force of rock dowels after support installation 


Concrete: 0.16m [ tegend 
Stag 


‘Support Element: plastic fiber shotcrete 


Figure 6.25 the safety factor of plastic fiber shotcrete 
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Hbeam(W): W150 x 18 [Togni 
. s 


Thrust N- MN 


Figure 6.26 the safety factor of steel rib 

6.7 Turning Bay (MAX. Buried Depth 150m) 

For turning bay, it is placed in rock class IIl (Its position is artificially controlled), 
according to LDP curve, the Characteristic curve presented as following, it is 


determined to choose 5m as the max. unsupported distance at support installation. 


Characteristic curve of Class III 


Total Displacement ui [m] 


0% 20% 40% 60% 80% 100% 
Stress Release Ratio Ai 


Figure 6.27 Characteristic curve 
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LDP Curve of Class III 
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Figure 6.28 Longitudinal Displacement Profile curve 
Some support system is considered: 
(1) The shotcrete strength is increased to 30MPa; 
The displacement without tunnel support, the displacement after support installation, 
plastic zone and axial force of rock dowels, the safety factor of plastic fiber shotcrete/ 


steel rib (including shear and moment) are shown as following. 


Total 
Displacement 
m 

0.00e+000 


4.00e-004 
8.00e-004 
1.20e-003 
1.60e-003 
2.00e-003 
2.40e-003 
2.80e-003 
3.20e-003 
3.60e-003 
4.00e-003 
4.40e-003 
4.80e-003 


3.28e-003 
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Total 
Displacement 
m 

00e+000 


-00e-004 
-00e-004 
-20e-003 
-60e-003 
-00e-003 


0. 

4 
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1 
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Z 
2.40e-003 
2.80e-003 
3.20e-003 
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4.00e-003 
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Figure 6.29 Total displacement without tunnel support 


Pa 2.20e-003 


2.80e-003 


Figure 6.30 Total displacement after support installation 
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Figure 6.31 plastic zone and axial force of rock dowels after support installation 
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Concrete: 0.08m [Tna 


Thrust N- MN 


Support Element: plastic fiber shoterete 


Figure 6.32 the safety factor of plastic fiber shotcrete 


6.8 Junction (MAX. Buried Depth 300m) 
The junction is inferred as rock class Il, according to LDP curve, the Characteristic 
curve presented as following, it is determined to choose 9m as the max. unsupported 


distance at support installation. 


Characteristic curve of Class II 


Total Displacement ui [m] 
o 
je) 
je) 
= 
UI 


0% 20% 40% 60% 80% 100% 
Stress Release Ratio Ai 


Figure 6.27 Characteristic curve 


Page 38 of 50 


DOOSAWN | UPPER TRISHULI-1 HEP (216MW) q tHg 


POWERCHINA 


Surrounding Rock Stability Calculation of Adit No.3 


LDP Curve of Class II 
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Figure 6.28 Longitudinal Displacement Profile curve 
The displacement without tunnel support, the displacement after support installation, 
plastic zone and axial force of rock dowels, the safety factor of plastic fiber shotcrete/ 


steel rib (including shear and moment) are shown as following. 


Total 
Displacement 


jä 2.73e-003 
0.00e+000 


3.00e-004 
6.00e-004 
9.00e-004 
1.20e-003 


1.50e-003 


1.80e-003 70e- 2.55e-003 


2.10e-003 
2.40e-003 
2.70e-003 2.85e-003 
3.00e-003 
3.30e-003 
3.60e-003 


Figure 6.29 Total displacement without tunnel support 
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Total 
Displacement 
m 

0.00e+000 


3.00e-004 
6.00e-004 
9.00e-004 
1.20e-003 
1.50e-003 
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™12.55e-003 


2.55e-003 


2.85e-003 


Figure 6.30 Total displacement after support installation 
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Figure 6.31 plastic zone and axial force of rock dowels after support installation 
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7 Wedge Analysis 


Concrete: 0.05m 


Support Element: plastic fiber shotcrete 


Figure 6.32 the safety factor of plastic fiber shotcrete 


7.1 Discontinuity Input Data 


According to the site investigation, there are three group joints may distribute in the 


adit tunnel, as is shown in Table 7.1. The stereographic projection of these joints in 


rock mass in shown in Figure 7.1~7.4. The Mechanical parameter of joints is shown in 


Table 7.2. 


Table hR CRS YA FB Ee RERANCH. .4 Occurrence of joints in adit tunnel 


Joint NO. Dip Dip Direction Tunnel Direction 
J1(foliation) 25° 333° 
J2 65° 160° 
N29°W 
J3 75° 120° 
J4 60° 65° 


Table #RRISCES BCA FRE PERCH. .5 Mechanical parameter of joints in the adit 


tunnel 
JCS C (0) 
JONIT TYPE JRC 
(MPa) (MPa) C) 
Foliation(J1) 6 70 0.20 35.0 
lithoclasts cemented 
8 50 0.15 28.8 
(J2, J3, J4) 
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Figure 4R CRY h AA H E PETRI CF. .1 Stereographic projection of the joins (J1/J2/J3) in the 


adit tunnel 
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Figure 4R h AA H E FERIE -2 Stereographic projection of the joins (J1/J2/J4) in the 


adit tunnel 
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Figure 4R CRY HCA H E FERAE -3 Stereographic projection of the joins (J1/J3/J4) in the 


adit tunnel 
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Figure 4R CRY h AA H E FERIE -o .4 Stereographic projection of the joins (J2/J3/J4) in the 


adit tunnel 


7.2 Wedge Stability Analysis 


According to the wedge stability analysis, potential wedges may occur in the 


excavation periods, which are shown in Table 7.3. 


Table 7.3 Instability Block Search Statistics of adit tunnel 
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Surrounding Rock Stability Calculation of Adit No.3 
Maximum 
Fissure Block a: Graphical 
block Safety | Development | Instability : 
S.N. set volume ; i presentation of 
height factor position mode 
number | (m°) block 
(m) 
D 
Lower left he 
1 J1fJ2/)3 0.291 1.63 0.824 f sliding l 
side wall 
Upper right N 4 
2 J1/32/)3 1.391 2.26 0.532 ; sliding 
side wall 
~ 
3 J1fJ2/)3 0.154 1.16 0.058 Crown fall 
————— 
| 
£> 
Ai 
Lower right May be yj | 
4 | Jifs2/s4 | 12.814 8.92 3.206 f A 
side wall stable A j 
Upper left E N N 
5 [312/34 | 38.11 12.26 | 0.361 ' sliding X 
side wall 
Aft ti hotcrete 50 4.156 MORE Q , 
er supporting, shotcrete 50mm ; ~ 
di £ stable 
M 
L 
6 J1/32/J4 | 0.087 2.51 0 Crown fall L 
Lower right Pen 4 
7 | 31/13/04 | 37.209 10.45 89.5 l sliding A í 
side wall /> 
CY, 
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LS | 
J 
May be 
After supporting, shotcrete 50mm 190 Le 
stable e 
k i 
Upper right e LA 
8 | J1/J3/J4 0.1 0.47 0.749 i sliding i 
side wall | 
&. 
Upper left a ~~ 
9 | 11/13/14 | 69.554 13.01 0.7 i sliding D 
side wall oN 
May be vt 
After supporting, shotcrete 50mm 2.99 
stable 
10 | J1/J3/J4 | 0.015 1.98 0 Crown fall j 
Lower left _ ko 
11 | J2/I3/J4 | 0.224 1.71 0 : sliding A 
side wall \ { 
A 
Upper right _ £ 
12 | J2/J3/34 | 0.306 1.99 0.195 ; sliding 
side wall / 
13 | 12/13/14 | 6.895 2.88 0.486 Crown fall y 
| 
, May be k 
After supporting, shotcrete 50mm 4.265 A 
stable 


It could be seen that, in all the possible joints combination except No.5/7/9/13, the 
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volume of the unstable block is about 0.1~1.3m?, and the height of the block is about 
0.4m~2.5m. In the tunnel, as for class | & 2, spot shotcrete can restrain the trend of 
sliding and fall. As for other classes, the pattern bolting can control the sliding and fall. 
During the excavation process, a reasonable layered excavation sequence shall be 
used and a reasonable excavation height to avoid large potential unstable blocks. 

As for the N0.5/7/9/13 joint combination, after spot and pattern supports, the 
minimum safety factor of potential wedges without supports in adit tunnel is bigger 


than the minimum safety factor of 2.0. 


8 Conclusions and Recommendations 


8.1 Conclusions 


As per the calculation fore adit tunnel above, it can be concluded as follows: 

(1) For rock class |, It can be seen that the maximum displacement is less than 0.1% 
of tunnel height (0.001x5600=5.6mm), the systematic support is not required. In the 
actual excavation process, spot shotcrete can be carried out according to the exposed 
rock. 

(2) For rock class Il, the tunnel roof occur plastic zone. In case of fall-block, the spot 
rock dowels many be required in the tunnel roof. And the thickness of shotcrete 
should be 50mm thick. 

(3) For rock class Ill, the initial proposed supports meets the requirement of tunnel 
stability, when enhance the strength of shotcrete, the thickness of shotcrete can 
remain 80mm for temporary tunnel. Doe to the adit NO.3 will be the access to 
headrace tunnel, the thickness of shotcrete increased to 100mm. 

(4) For rock class IV, while the maximum buried depth is no greater than 100m, the 
initial proposed supports meets the requirement of tunnel stability; while the 
maximum buried depth is greater than 100m and less than 150m, the thickness of 


shotcrete should be improved to 120mm and the strength of shotcrete should be 
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improved to 30MPa and the rock dowel spacing is decrease to 1m; while the 
maximum buried depth is greater than 150m and less than 200m, more measures 
should be provided, e.g.: The shotcrete thickness is increased to 160mm, the 
shotcrete strength is increased to 30MPa, the spacing of Steel support and rock dowel 
is decreased to 1m. 

(5) For rock class IV, only the buried depth below 200m is considered, and special 
study is conducted when the buried depth exceeds 200m. 

(6) At the turning bay and junction transition section, the strength of shotcrete 
enhance to 30MPa. 

(7) For wedge analysis, after shotcrete, the minimum safety factor of potential 
wedges without supports in adit tunnel is bigger than the minimum safety factor. 
However, the drainage system in construction periods is extremely important in the 
stability calculation of wedges, it should be paid more attention in the excavation 


process. 


8.2 Recommendations 


According to the calculation above, the recommended supports for each rock type are 


summarized in the Table 8-1. 
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Table RIA PRA HEERE. -6 Recommended supports for each rock type 


Initial proposed supports 


Item Unit 
| j Junction Il | Turning IV 
Max. buried depth m 300 |300 300 300 150 200 |150 |100 
Shotcrete thickness mm Spot | 50 50 100 80 160 | 120 | 100 
Shotcrete strength MPa 25 25 30 25 30 30 | 30 | 25 
Rock dowel length m 3.0 3.0 3 3.0 3 3.5 | 3.0 |3.0 


Rock dowel spacing within 


m spot | spot 2.0 2.0 2.0 1.0 1.0 1.5 
each row 


Rock dowel row spacing m spot | spot 2.0 2.0 2.0 1.0 | 1.0 11.5 


Steel support/Lattice girder 
spacing 


Stability Yes Yes Yes Yes Yes Yes | Yes | Yes 
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